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Abstract 

The Bird and Nanbu systems are particle systems used to approximate 
the solution of the mollified Boltzmann equation. In particular, they have 
the propagation of chaos property. Following |GM94| . we use coupling 
techniques and results on branching processes to write an expansion of 
the error in the propagation of chaos in terms of the number of particles, 
for slightly more general systems than the ones cited above. This result 
leads to the proof of the a.s convergence and the central-limit theorem 
for these systems. In particular, we have a central-limit theorem for the 
empirical measure of the system at a fixed time t under far less stringent 
assumptions then in [M_ei9_8j. 
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1 Introduction 

In a recent work ( |DPR09aj l. we showed a expansion of the propagation of 
chaos for a Feynman-Kac particle system. This particle system approximates a 
particular Feynman-Kac measure, in the sense that the empirical measure asso- 
ciated to the system converges to the Feynman-Kac measure when the number 
of particles N goes to oo. What is called propagation of chaos is the prop- 
erty of the particle system that q particles, amongst the total of N particles, 
looked upon at a fixed time, are asymptotically independent when N +oo 
{q is fixed) and their law is converging to the Feynman-Kac law. In [DPROQa] . 
we wrote an expansion in powers of N of the difference between the law of q 
independent particles, each of them of the Feynman-Kac law, and the law of q 
particles coming from the particle system. One can also call this expansion a 
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functional representation like in |DPR09a] : in the present paper, we call it an 
expansion of the error in the propagation of chaos. In the setting of |DPR09a] . 
the time is discrete. We showed there how to use this kind of expansion to 
derive a.s. convergence results (p. 824). In ( |DPR09b) ). we extend the result 
of |DPRQ9aj to the case where the time is continuous, still in the Feynman-Kac 
framework, and we show central-limit theorems for {/-statistics of these systems 
of particles. The proof of the central-limit theorems for {/-statistics relies only 
on the exploitation of the expansion described above. 

We wish here to establish a similar expansion for a family of particles systems 
including Bird and Nanbu systems (that is, the equivalent of Theorem 1.6 and 
Corollary 1.8 of |DPR09bj ). Bird and Nanbu systems are used to approximate 
the solution of the mollified Boltzmann equation. We refer mainly to |GM97j 
and take into account models described in (2.5), (2.6) of [GM97| (a similar 
description can be found in |GM99) . Section 3). An other reference paper on 
the subject is |GM94) . The two main points of interest of this paper are: it 
provides a sequel to the estimates on propagation of chaos of |GM97j . |GM99j 
and it allows to apply the results of |DPR09a] . |DPR09bj to Bird and Nanbu 
systems. In particular, 

• we get a.s. sure convergence (Th. 15. ip (completely new, based on our 
knowledge) 

• and we get a central-limit theorem for the empirical measure of the system 
(Th. EjI) with far less stringent assumptions then in |Mel98j Th. 4.2, 4.3. 

The proofs leading to the development in the propagation of chaos are radically 
different from those in |DPR09b] and this is why we decided to write them in a 
different paper. 

In Section O we will recall the definitions of Bird and Nanbu models, as can 
be found in |GM97) and will give an equivalent definition, useful to our purposes. 
In Section [3l we will state and prove our main theorem about the expansion of 
the error in propagation of chaos (Theorem 13. ip . The proof relies on estimates 
on population growth found in |AN72) and on coupling ideas. In Section [H we 
prove in what is called a Wick-type formula in iDPR09aj (see (3.6) p. 807 in 
|DPR09aj and |DPR09b j. p. 15 and PropositiongSl, this formula (Corollary[431 
and Proposition 14.31 are used in Section [5] to prove an a.s. convergence theorem 
for the empirical measure (Th. 15. ip and central-limit theorem for the emprirical 
measure (Th. [OJ- 

2 Definition of the model 
2.1 Bird and Nanbu models 

In all the following, we deal with particles evolving in E := M'*. We set the 
mappings e, : h £ R'^ ^ e,{h) = (0, . . . , 0, /i, 0, . . . , 0) € K''^^ {h at the i-th 
rank) (1 < i < N). We have a kernel ^{v, w, dh, dk) on R^'' which is symmetrical 
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(that is , w, dh, dk) = 'fi{w, v, dk, dh)). We set iJ,(v, w, dh) to be the marginal 
'jj,{v,w,dh X M.'^). Our assumptions are the same as in |GM97] : 

(HI) We are given a Markov generator L on R"* acting on a sufficiently large 
domain V{L) of L°°(R'^). (See |GM97| p. 119 for a discussion on V{L)). 

(H2) We suppose sup^ „ /i(x, a, M'' x M'') < A < oo. 

In Nanbu and Bird systems, he kernel /t and the generator L have specific 
features coming from physical considerations, the coordinates in R'' represent 
the position and speed of molecules but these considerations have no effect on 
our proof. This is why we claim to have a proof for systems more general than 
Bird and Nanbu systems. 

The Nanbu and Bird systems are defined in (2.5) and (2.6) of [GM97J, by the 
mean of integrals over Poisson processes. We give here an equivalent definition. 

Definition 2.1. The particle system described in ]GM97^ is denoted by 

{Zt )t>0 — (Zi ' )t>OS<i<N ■ 

It is a process of N particles in W^' and can be summarized by the following. 

1. Particles (.^^'*)i<i<Ar in R"^ are drawn i.i.d. at time according to a law 

2. Between jump times, the particles evolve independently from each other 
according to L. 

3. We have a collection (-/Vi.j)i<i<j<jv of independent Poisson processes of 
parameter A/ {N — 1). For i > j, we set Nij — Nj^i. If Ni_j has a jump 
at time t, we say there is an interaction between particles i and j and we 
take a uniform variable U on [0,1], in dependant of all the other variables, 

if U < — - then the system undergoes a jump: 

• In the Bird system: = + ei{H) + ej{K) with 

/I(^fJ^zfJ^R2^) 

(independently of all the other variables). 

• If we replace p. by fJ.'{z, a, dh, dg) — fi{z, a, dh) (g) 5Q{dg) + 6o{dh) (g) 
fi{a,z,dg) in f2.1]) . we obtain the Nanbu system (cf. Remark 2.6, p. 
120, ]GM97j ) 

Theorem 3.1 of |GM97j implies that there is propagation of chaos for this 
system. This theorem says (Vg,t): 

||£(zf^...,zf-)-/:(zf>|i..<2,(,-i)^^ + ^'^' 
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and 

|l£(zf'Vi^*llTy<6^^ , 
where (Pt) is solution of (with Pq fixed) 

V<^ e viL), 
dt{^,Pt)~{L^,Pt) 

-(0(z + h) — <j){z) + (t>{a + k) — 4){a))fi{z, a, dh, dk), Pt{dz)Pt{da) 

We can deduce propagation of chaos from the previous results, that is Vt, VF 
bounded measurable, 

\ . . . , Zf-)(F) - PriF)\ < (2,(, - 1)^^^ + 6^ 

In Theorem 13. H we will go further than the above bound by writing an expan- 
sion of the left hand side term above in powers of N. We will use techniques 
introduced in |GM97j . The main point is that we look at the processes backward 
in time. 

2.2 Backward point of view 

From now on, we will work with a fixed time horizon T > 0. For any j G N*, we 

set [j] = {1, . . . , j}. For A > 0, we call £{X) the exponential law of parameter A. 
We fix g e N*. We start in s = with = {i}, \/i e [q]. We set Vi, 
^ ^C*. For 1 < i < j < A^, we define processes {Ni'^ = iV|'*)s>o and 

for i e [q], we define (C|)s>o, (^^)s>o (respectively in N, V{N), N*) by the 

following. For all s G [0,T], we set 

The processes {N^'^), (C*), (K^) are piecewise constant and make jumps. We 
define the jump times recursively by (taking {Uk)i<i<q.i<k, {Vk)i<i<q,i<k i-i-d. 
^ 8{l))To = and (always with the convention inf = +00 and (■■•)+ standing 
for the nonnegative part). 



Tk = inf i^Tk-i < s <T : J _ ^ "^^^ du > Uk 



Tk = inf(r^,r/'^ 

In Tk-. 
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• If Tfc = T^., we take r(fc) uniformly in C^^,- U • • • U C!^^_ and j{k) uni- 
formly in [N]\{C^^_ U • • • U C^^_). For any i such that that r(fc) € C^^_, 
we then perform the jumps: C\ = C\_ U and iV''('=)'J('=) (Tfc) = 
jVr(fe),i(fc)(T^_) + i_ 

Notice that the (•■•)+ the definition of above forbids to be in the 
situation where we would be looking for j{k) in 0. 

• If Tfc = r^', we take r{k) uniformly in C^^_U- • •UC|.^_ and j{k) uniformly 
in C^^_U- • •UC|.^_\{r(A:)}. Foranyi,z' such that r(A;) e C^^_ andj(fc) e 
C^^_, we perform the jumps: C^^ = C|.^_ U {j(fc)}, C'^ = Ci^^_ U {r{k)} 
and 7V('=)j('=)(Tfe) = iV'^WjW(rfc-) + 1. 

For all i, i, we set ifj — #6*4. This whole construction is analogous to the 
construction of the interaction graph found in |GM97j . p. 122. 

We now define an auxihary process (Zi^)o<s<T = {Z^''^)o<s<T.i<i<N of N 
particles in M''. 

Definition 2.2. The interaction times of the {Z^'^)i<:s<T,i<i<N are {T — 
Tk,k > l,Tk <T}. (We say that the times (T^) are defined backward in time.) 

• Zq , . . . , Zq are i.i.d. ^ Fq 

• Between the times {T — Tk)k>i, the Z^'^'s evolve independently from each 
other according to the Markov generator L. 

• At a jump time T — Tk where is a jump time of N^'^ , (Z^) undergoes 
an interaction having the same law as in Definition \2.1\ ^). 

Definition 2.3. For all t>Q, we set 

Lt = MkeN:Tk<t,Tk^T;!} . 

We call this quantity the number of loops on [0, t]. 

Example 2.4. Take q — 2. Suppose for example, that the only jumps of the 
7V''^ 's occurring in [0, T] are 

AN^'^{2T/3) = 1 , AiV^'2(T/3) = 1 

then 

. for s e [0, T/3[, = 2, = 0, K] = = 1, 

• for s e [T/3, 2T/3[, K, = 2, = 1, K] = = 2, 

• forse [2T/3, T], K,^ 3, L, = 1, ^ 2, ^ 3 . 

We have to keep in mind the following lemma throughout the whole paper. 
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Lemma 2.5. 1. //Yi, ■ ■ ■ ,Yk are independent and respectively of law £{Xi] 
. . . ,£{Xk) (Xi , Xk >0, k€W) then 

mf(Yi,...,Yk)^£iXi + --- + Xk) , 

Ai 



p(ri = inf(yi,...,n)) = 



Ai + • • • + Afe 

2. For Y - £{X) (X>0), C{Y - t\Y >t)^ £{X) (for any t>0). 

3. IfYi,Y2,... i.i.d. - £{X) (X > 0), k e W, Ui,U2,...,Uk is the order 
statistics of k uniform random variables on [0, T] then 

L{{Yi,Yi+Y2,...,Yi + --- + Yk) 
\Yi^--- + Yk <T <Yi + --- + Yk+i) 
= Ci{Uu...,Uk)) 

We then have: 

Lemma 2.6. For all T>0, {Z^-\ . . . (Z^'\ . . .,Z^''). 

The system {Z^)q<:s<y is of use in Section|4]but is also useful to understand 
the next auxiliary process, which we use in Section[3l We now define, for a fixed 
time horizon T > 0, the auxiliary process (Zj^)o<5<T — (^o<!,<T)i<i<w,o<s<T- 

We start in s = 0, with Cq = [q], Kq ~ q, Lq ^ 0. We define processes 
(Cs)s>o, (i^)s>o, {Ls)s>o (respectively in ^'(N), N*, N*). These processes are 
piecewise constant and make jumps. We take {Uk)k>o to be i.i.d. ^ £{l) and 
(^fe)fc>o to be i.i.d. ^U{[0,1]) (these variables are independent of all the other 
variables). We define recursively the jump times (Tfe)fc>o by Tq = and 

Tk+i = mf . > Tk : AK^du > Uk 



InTfc: 



• If Ak < then we perform the following jump: K^p = 

Kf^_ + 1, Lf^ = Lf^_, we choose i{k) uniformly in [N]\Cf^_ and Cf^ = 
Cjr_ U {i{k)}. We choose j{k) uniformly in Cf^_. 

• If ^fc > ,.r^^^~,^"''i> then we perform the following jump: Kf = 

Kf^_, Lf^ = Lf^_ + 1, Cf^ — Cf^_. We choose i{k) ^ j{k) uniformly in 

-I fc — 

Definition 2.7. The interaction times of the (^i^'*)i<s<T,i<i<A' are {T — 
Tk,k > l,Tk < T} . (This is why we say that the times (Tk) are defined backward 
in time.) 
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• Zq , . . . ,Zq are i.i.d. ^ Fo 



• Between the interaction times {Tk)k>i, the Z^'^ 's evolve independently 
from each other according to the Markov generator L. 

• At an interaction time T — Tk, {Z'^) undergoes an interaction having the 
same law as in Definition ] 2. 11 with i,j replaced by i{k), j{k). 

Keeping in mind Lemma [275} we get: 
Lemma 2.8. For all T>0, 

(Kt, Lt)o<t<T [Kt, Lt)Q<t<T ■ 

3 Expansion of the propagation of chaos 

We define for any N,q e W , q < N: 

Nl 



{q,N) = {a : [q] ^ [iV], a injective } , (TV), = #{q,N) = 

l<i<N 



"i {N-q)\ 

Let us set 

N 



ae{q,N) 

For any function F : M, we call (r/f a [/-statistics. Notice that 

for all function F, 

E{F(Z^'\ . . . X'")) = E((^f . (3.1) 

We define 

-Psym J ■ ■ ■ I 2:') — — ^ ^ F(^X ^ \ . . . , X ^''^ ) , 

where the sum is taken over the set Sq of the permutations of [q] . Notice 
that 

(r?f )°'(F) = (r;f )°«(F,,^) . 
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Theorem 3.1. Set a — e For all q > 1, for any bounded measurable F, 
VT > 0, yia > 1, 

^iiVrfHF)) = Pt''{F)P{Lt = 0) (3.2) 



E 

l</<io 



1 -^N,lo + l 



where the A^^, A^j."^^ are nonnegative measures uniformly bounded in N 
defined by, for any bounded measurable F, 

A^r^iF) = E{FiZ^-\ . . . , Z^-^)\Lt = ^ l)iN - 1)' 

A^JriF) = . . . , Z^^«)|iT > OlP(iT > 1){N - 1)' . 

K^e further have the following bounds (VF £ C^iW') ) 

sMAl:^{F)X:kF)) 

qil^afl^-^ (2/ + 1)! 

- n (l-(l-«)l/9)2i+2 ^ ll^ll- 

(1 - (1 - 7V>i ^ ^ 

Let us define P^.^C^) = "^{{tIt )®''(^))- Using the terminology of |DPRn9a,) . 
p. 782, we cannot say that the sequence of measure i^T,q)N>i is differentiable 
up to any order because the A^^ appearing in the development depend on N . 

For a specific class of functions F, we will compute what is the order of 
^^t'^^^ (F) (Proposition gJl Corollary g^]). This last result, together with the 
boundedness of the A^^, is what we need on the proofs of Section [H 

Proof. According to |GM97) (section 3.4, p. 124) or, equivalently |GM94) (sec- 
tion 5), 

E(F(Z^'\ . . . , Z^'')|Lt = 0) = Pt{F) . 
We have by Lemma [TBI V/q: 



E(F(Z^'\...,Z^'')) = E(F(Z^'\...,Z^'')|Lt-0)P(Ft = 0) 



la 



1=1 

+E(F(Z^'\...,Z^''«)|Lt>/o + 1) 
xP(Lt> ^0 + 1)- 

It is sufficient for the proof of l|3.2p to show that P(Ft > is of order < 

V/ e N*. 
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We define piecewise constant processes {K[.)s>o, iL'g)s>a (in N) such that 
K'q — q, L'q ~ 0. Their jump times are (r^.)fe>o defined recursively by Tq = 
and 

f'^infJs: / AK'du>uA • 



InT^: 

• If^fc < ]v3i , then we perform the following jump: if-, — K)~, +1, 

K n- 

• If ^fe > , then we perform the following jump: Ki-, — K~ +1, 

L' = L' + 1. 

Notice that we use there the same variables UkS, and A^'s coming from the 
definition of {Z^). We have: 

w, 2{N-kt{io))+ + kt{Lo)-l 
as Vt, uj, 2{N -\) *■ 

then Vi, ^t(w) < k'^[uJ) , 

2(iV-^«(a;))+ (TV - j?;(c.)) + 
as Vr, uj, > — — ; 

2{N~Kt{Lo))+ + Kt{Lo)-l - N-1 

then \/t, UJ, Lt{uj) < L'^{uj) , 
yt, Lu, L[{uj) < k[{uj) . 

The process {K'^)s>q is equal in law to the sum of q in dependant Yule 
processes Ys^\ . . . , y}'^^ and its law is thus independent of N (see |AN72j . p. 
102-109, p. 109 for the law of the Yule process). We have V{yP = k) = 
e-'^{l- e-^^ f^^ and so: 

v{k[ ^k) = P(y/^^ + • • • + r/'^' = fc) 

< ^p(r«>rfc/gl) 

i=l 

< q{l - e-'^fli-^ . (3.4) 

Notice that Vi e [0, T\ 

P(AL; = 1|AK,' = 1, (x;)t>o) = 1 - < ^ . (3.5) 
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We decompose 

F{Lt = < ViLr > I) 

(by Lem. EH = P(Lt > 

< p(£t > 



< V{K'j.>lVW^\)+¥{L'j,>l,K'j,<[^/W^\) 

and we compute 



P(L^ > l,K!j, < [VA^- IJ) 

= ^(^T = r\K^ = /c)P(^^ = fc) 

/<r<L\/lV-^J r<fc<L\/lV-^J 

i<T<Lyiv-^j r<fc<Lyiv^j 

;<fc<|_Viv-^j '<''<fc 
As, for fc < [ViV - IJ , < 1, we get 

i<fe<L\/lv^j 



^ /!(jV_iy Efe(fe + l)(fe + 2)---(fe + 20 

x((l-a)i/9)fe-i(i_Q,)i/<7-i 
(2/ + 1)! 



< 



;!(7V-1)' (1 - (1 - a)!/-?) 



2;+2 



We also have 



< 



g(l_a)Lv^J/g-i 

(1- (l-a)i/9) 

q 

(1- (1 -a)i/9) 

x»p(,i-„.)iv^i/»-.(»-iy)^ 



□ 
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4 Wick formula 



We now define an auxiliary system {Zl)o<t<T,i>i with an infinite number of 
particles. We start in s = with (5^ = {i}, Vi e [q]. We set Wi, K^^^ #(7^. For 
1_< j < j < A^, we define processes (iV]'-' = iV|'')5>o,i jeN* , (C|)s>o,i<j<g, 
(i4r*)s>o,i<i<g (respectively in N, P(N),N) by the following. The processes 
(iV''^), (C*), {K^) are piecewise constant. We set 

We define the jump times recursively by (taking {Uk)i<i<q,i<k i-i.d. ^ ^(1)), 
To = and 

ij, = i.f {f._. <.<T:f_ KK.- A^..(;;;^-..). ,„ > g, j 

ffe = inf(r^,inf{Ti >ffc_i}) 

(recall that the process [Kt) and the T^'s are defined in Subsection l2.2p . Notice 
that {Tfc, /c > 0} C {Tfc, A: > 0}. In ffe: 

• lifki {T{', I > 1}, we take r{k) uniformly in C^^_ U • • • U and J(fc) 
uniformly in N*\(C'^^_ U ■■■UC^^_). For any i such that r{k) G C'^^_, 
we perform the jumps: — C'^p _ LI {j{k)}, = -^t - + 1 ^^d 

• If Tfe e {T/', / > 1}, we take r{k) uniformly in C^^_ U • • • U C^^_ and J(fc) 
uniformly in C^^_U- • ■UC!^^_\{r{k)} . For any z, i' such that r(fc) e C't^-j 
j{k) G C^j^_, we perform the jumps: = C^^_U{j{k)}, K^^ = K^^_ + 
1, = U {r(k)}, = + 1 and ivj''^'^'*''^ = TvJ*!^'^'^''^ + 1 

Notice that Vi e [g], the processes iiT* and #C" are not equal. The following 
lemma is a consequence of Lemma 12.51 

Lemma 4.1. 1. The process {Ks)s>o is piecewise constant, has jumps of size 
1 and satisfies \fO < s < t 

And so it has the same law as {K[.)s>o. 

2. For all t, Kt>Kt, a.s. 

3. //Ti = Ti, . . . , Tfe = ffe then Kt, = Kt,. 
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Definition 4.2. The interaction times of the are {T — Tfc,/s > 1} (we say 
they are defined backward in time). 

• The (Zl) are i.i.d. - Pq. 

• Between the jump times, the evolve independently from each other ac- 
cording to the Markov generator L. 

• At a jump time T — Tk, (Z) undergo a jump like in Definition \2.1i {^j, 
with i,i replaced by r{k),j{k). 

We define (Vt > 0) the event G and the trajectories JCt, JCt by 

G = {Vfc > 1 such that fk <T, f^^ T^} 
JCt = {{K)l<s<t,^e[q]} 
ICt = {iKl),<,<t,te[q]} 

and, for q even, we set {Tk^ < Tk^ <•••} = {T^ : fc > 1, = T^'}, {T^^ < 
Tj, <...} = m:fc>l,3?,ffe=T/'}, 

A = {#{fc> 1 :Tfe =r^',Tfc <T} = g/2} 

n{r(fci)eC^,^_,j(fci)eC|,^_,... 

A = mkJ^l:3l,n~Tl',fk<T} = ql2} 

n{F(fci)eCi.^^_,j(fci)eC|^^_,... 

...,^(^,/2)eC|:i _J(fc,/2)eCi J 

(recall the i{k)^s, r(fc)'s are defined in Subsection I2.2p . For q S N*, we say that 
F : (M'')'? — > R is symmetric if for all a in the set of permutation of [q] (denoted 
by Sq), Vxi, . . . e M'', F(a;^(i), . . . = F(xi, . . .,Xq). We define 

BQ^"^{q) ~ |f : {W^y ^ M, F measurable, symmetric, bounded, 

F{xi,...,Xq)PT{dXq)^0} . 

Proposition 4.3. For F G B^"^'{q), we have: 

• fork< f , 

• for q even. 



for k < f , A^^^'iF) = 0, 



iV'/2E((,7^)®''(F)) _^_^E(E(i^(Zi,...,Z^))|/C„A) 



w^+oo 2«/2(g/2)! 



n f KKf~'Kfds) . {A.l) 



X 

\<i<ql2 
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Proof. In the following computations, we set R — sup{fc : Tk < T}. We have 
E(G|/Ct) = E(lj,,^^ J/Ct) 



^Ti=Ti ■ ■ ■ ^Tr=Tr' 

xE(i^^^^^|/Ct,Ti = ri,...,rj^_i = Tfl_i)) . 

On the event {Ti — Ti, . . . , Tr-i = Tfl_i}, we have Kt — Kt, yt : Tr-i < t < 
Tp>. And so we have: 



I ^T,Ti =Ti, . . . , Tfl_i — Tfl_i 



>PI / A^rf^<[/, 



2 



= exp ( — {Tr - Tk-i) 
So, by recurrence, 

E(G|/Ct) > exp [-J^^T^ ■ (4-2) 

When 3fc, r, s : A^'^''" jumps at s, fc e C*, r e and s € [0, T], we say there 
is a loop at s between C" and . If i — j, we say that has a loop with 
itself at s. We can define in the same way loops between C*, Notice that 
[3i, C has two loops on [0, T] and Lt < q/2] => [ 3j, has no loop on [0, T]]. 
Conditionally on {C^ has no loop on [0,T]} (i.e. has no loop with any other 

and has no loop with itself), Zj,'^ is independent of . . . , Zj, '"^ and 

has the law Pt (by Theorem 13.11 applied for g = 1 and Lemma [TS]) . Suppose 
k < q/2, as G Sq^™((j'), we then get 



E{F{Z^'\...,Z^''')\LT^k) 

= E(E(F(Z^'\...,Z^'«)|Ci has no loop on [Q,T], Lt ^ k)\LT ^ k) 
= (4.3) 

We suppose by now that q is even. Because of Theorem 13.1) we then get: 

Ar9/2]E((^A')0?(i^)) ^ N'^/^E{F(Zt\...,Zt''^)) 
= Ar9/2E(F(Z^^\...,Z^^'?)) 

which has same limit as Af9/2E(F(Z^'\ . . . , Z^''')lz,^^^/2) when A^ ^ +oo. As 
2<!/^q/2y. number of ways of partitioning [q] into q/2 couples, we get 
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Noticing that £((Z™'\ . . . , Zp'^)\G) = C{{Z^'\ . . . , we get: 



2«/2(q/2) 



-E{F{Z^'\...,Z^''')1^1g 

Notice that, knowing JCt, for i j, the number of foops between and C-' on 
[0, t] is a non-homogeneous Poisson process of intensity ( ^wlT' 1 ■ Notice 

V / 0<t<T 

also that, knowing /Ct, for i ^ j, the number of loops between C* and on 
[0,t] is a non-homogeneous Poisson process of intensity ^ ^ . So we 

have: 

P( no loop between C" and C^|/Ct) = exp I - y ~/v^~l^^'^ ) ^' ' 

~ - ~ ( r AK'-K^ \ 

P( no loop between and C-' |/Ct) = exp I - y ~/v^l^^* ) • 

We set 

B = {at least one loop between and C^} fl . . . 

• • • n {at least one loop between C'~^ and C'} . 

We have 

F{A\ICt) = P{B\ICt) ~ P{B\A\ICt) ■ (4.4) 

and 

P{B\ICt) = n (l-a(2j-l,2j)) 

l<J<9/2 

^ ^ l<J<9/2 

(byLem.SU) < [] (TAK^'^'K^t') ■ (4-5) 
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So, we have: 



N'^/'E{\{F{Z^-\...,Z^''')U1g^\) 

< m^''\\F\\ooE{E{lAlG^\JCT)) 

= 7V«/2||f|UE(E(1^|/Ct)E(1gc|/Ct)) 
(using {421), 63), dlSl)) 

(using Lemma |4?T|) 



(A^- 1)9/2 



X 1 — cxp 



K{KtYT 
N -I 



For a fixed 



^ ^ l<j<q/2 

— > (4.6) 



and 



i<j<g/2 

which is of finite expectation by l|3.4p and Lemma ICT So, by dominated con- 
vergence: 

m'^n\{F{z^-\...,z^^'')\A\GA) 0. 

We can show in the same way : 

iV«/2E(|(^^(4,...,Z^)l^lGe|)^ 0. 
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We have: 

7V9/2e(F(Z^,...,4)1^) = E(E(F(4,...,Z«)|^T,A))Ar9/2p(I|^T)) 



We set 



B = {at least one loop between and C^} D . . . 

• • • n {at least one loop between C^~^ and C^} . 



We decompose 



V{A\Ct) = P(-B|/Ct) - P{B\A\ICt) • 



We have 



iV«/2p(S|/CT) = Ni/^ H (l-5;(2j-l,2j)) 

l<i<9/2 



N- 



+00 ^ ^ /r, 

1 /o 'J 



l<i<q/2 



and 



P( at least two loops between and C^|/Ct) < (1 - "(1,2))^ . 
So we get, computing very roughly: 



¥{B\A\ICt) < E E 

l<r<q/2 l<ii,...v<g/2 



n (l-5(2j-l,2j)f 

je{ii,...,ir} 



X [] (l-5(2j-l,2j)) 

ji^{n,...,ir} 

2 



l<r<g/2 l<j<g/2 



[N - 1)9/2+1 



< 



29/2 



(iV - 1)9/2+1 

l<J<9/2 '- 



SO 7V9/2P(B\A|/Ct) ^ 0. And so: 

JV— >+oo 



JV— >+oo , /a 



l<i<g/2" 
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and again by dominated convergence, we get the result. 



We denote by Tq the set of partitions of [q\ into pairs. 



□ 



Corollary 4.4. [Wick formula] For F e B^^'^iq) of the form (/i (g) • • • (8) fq)sym 
and q even, 

with ^ 

/,) - E(/,(Zi)/,(Z|)|I, I,) X / KKlKlds . 



5 Convergence theorems 

All the theorems of this section are valid under assumptions , (121) • 

5.1 Almost sure convergence 

Theorem 5.1. For any measurable hounded f , T > 0, 

^t(/)„^ Prif)- 

N—*+oo 

Proof. We recall the notations of [DPROQci] . For any empirical measure m{x) = 
(based on N points x^,x'^, . . . , x^), any q, 

'^(^)^'' ]^ ^(:cMi),...,x»(<!)) , 

ae[Af][9l 

where [A^]M ^ a : [q] ^ [N]. Notice that for any F 

We define, VI < p < q, [qp := {a G #Im(a) = p} and (Vfc < q) 

(the s(., .) are the Stirling numbers of the first kind) and Vi^(of q variables), 
i?b(i^)(xi,...,x«) = F(x''W,...,a;^(«)) , 
Do^L,{F)= Y s{p,q-k)j^ Y F{x<'\--- ,x'^^'^^) . 
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We then have, by Corollary 2.3 p. 789 of |DPR09a) . for any empirical mea- 
sure m(x) (based on N points), any F of g variables, 

\o<k<q J 

Suppose F e BQ^"\q), 

E{{7j^r^F))^ E ]^ E <p^<i-k) E ^Hv^r'iDaiF))) . 

0<k<q q~k<p<q oSMI?' 

Notice that for 0<k< q, q~k<p<q, ae [q]^p\ #{i G [q],#a{{i}) = 
1} > q — 2k. Suppose by now that q is even. If fc < q/2, using Theorem 13.21 and 
reasoning in the same way as we did to prove (|4.3p . we get (for some constant 
C which will change from line to line): 

E{irj^r^DaiF)))<j^. 

And so 

E((^^)^^(F))<-^. 

iV 2 

Suppose we take a bounded function f : E ^ R. We set f = f — Pt{I)- We 
then have (with the notation /®''(x\ . . . , x') := /(a;^) x • • • x /(a;')) 

m{iT{f)-PT{m = niv^im 

= E((r7^)«'' (/««)) 
C 

< — ■ 

Ni 

Provided we take q big enough, we can apply Borel-Cantelli Lemma to finish 
the proof. 

□ 

5.2 Central-limit theorem 

Theorem 5.2. For aUfi,...Jqe 6o^'"(l), VT > 0, 

7V^/2(^^f(/i),...,^^(/,))^^ Af{0,K), 
with K{i,j) = Mhfj) + v^ih. fj). 
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Proof. For any u\, . . . ,Uq, we have 
E ( exp [Nr]^ (log ( 1 



Af / , / 1 , H 1- iUqfq 



= ^ (^""P ^^P^A^'-'/^T ii^Ulh +■■■+ iUqfqf) 

~ E(exp(x/]V(mi»7^^ (/!) + ••• + (/,)))) 

JV— >+(x) 

xexp('-ii3T((wi/i + ---+V9)'))) ■ 



We also have 
E (exp (^7Vr/|f I log 



^ _^ Wl/l + . ..iUqfq 



\/]V 



\0<fe<JV l<ji,-,jfc<« 



l<ii<---<ifc<A'' 



M E 



0<k<N l<j 



J2 ^'^^j.■■■Uj,ll{v?r''{fn^■■■^fj.)] 

i,-jfc<<? " / 

.-oo E (-1)'=^^ E ^^^E n 

fe>0,fe even l<jivJis<9 ^fceXfe {a,6}e/fe 



2'=/2(fc/2)! 

K>0,fc even ^ ' ^ 

fc/2 

E 2^74/2)! ( E (/,„/,,)) 

fe>0,fc even \ I / \\<j\,j'2<q 



= ^''P ~^ E "ji"j2E(l^#(/ji, /jj) j 

\ l<jl,j2<9 / 



□ 
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